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(54) Method of predicting and optimizing fuel production 



(57) A process for the prediction and the optimiza- 
tion of the output of a plant producing products from in- 
coming materials. The incoming materials are classified 
according to various physical characteristics and costs. 
The desired or ordered products arc also classified ac- 
cording to price and physical requirements. The incom- 
ing materials information and the product information is 



entered into a database which is accessible by a com- 
puting device. The computing device then calculates the 
optimal production process by calculating a plurality of 
production cycles and selecting the cycle with the opti- 
mum profitability. The computing device is programmed 
with non-linear equations derived from a regressive 
analysis of data collected from samples of incoming ma- 
terials and products. 
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Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

[0001] The present invention relates generally to a process and method for the prediction of the properties of and 
the optimization of a plant's output of products from a source or sources of raw material. 

[0002] More specifically, the present invention relates to a process and method for increasing the predictability and 
10 profitability of operations where a series of raw materials are combined and processed into intermediate or final products 
by optimizing the cost structure of the raw materials, and the output of final or intermediate products to result in the 
lowest cost materials input and highest value production output. 

[0003] The present invention also relates to the optimization of refining processes and petroleum blending operation 
to result in the highest value production output from available fuel stocks. 
i5 [0004] The present invention also relates to the accurate prediction of final properties of a blended fuel utilizing non- 
linear optimization and property prediction. 

Description of the Prior Art 

20 Fuel Additive Industry 

[0005] It is well known in the petroleum and other industries to maximize profitability by blending fuel stocks with 
appropriate additives to increase their commercial value. Ethyl Corporation, of Richmond Virginia supplies many fuel 
additives, including a diesel cetane improver known commercially as DII-3 rM which is used to raise the cetane level of 

25 a diesel fuel stock and thereby make otherwise lower-valued fuel stocks into valuable commercial fuels. MMT (man- 
ganese methylcyclopentadienyl manganese tricarbonyl) is a fuel additive, also manufactured by Ethyl Corporation, of 
Richmond Virginia, that provides octane enhancement while reducing the amount of crude oil necessary to produce 
gasoline. Ethyl Corporation additionally manufactures the HiTEC® 4700 Series of antioxidants, including hindered 
phenolics and diphenylamines to JP8+100 jet fuel additive which have known effects on the fuel. 

30 [0006] Performance fuels for varied applications and engine requirements are known for controlling combustion 
chamber and intake valve deposits, cleaning port fuel injectors and carburetors, protecting against wear and oxidation, 
improving lubricity and emissions performance, and ensuring storage stability and cold weather flow. 
[0007] Fuel detergents, dispersants, corrosion inhibitors, stabilizers, oxidation preventers, and performance additives 
are known to increase desirable properties of fuels. 

35 [0008] It is known that mixtures of fuels and additives can increase and decrease desired properties in a resultant 
fuel blend. 



SUMMARY OF THE INVENTION 

40 [0009] The present invention contemplates supplying to an end user a tool for the prediction of and the consequent 
optimization of production from a plant which creates products from available raw materials. 
[0010] The tool includes an update-able database of incoming raw material properties, preferably including physical 
properties as well as economic properties, e.g. the cost of those raw materials. The tool also includes a database of 
blendstock properties, product specifications; processes, and the market price for the final and/or end product. A user 

45 can input the available raw material stocks, including their costs and available volumes, and input the standing pro- 
duction orders, including price and minimum acceptable values, and the computing device will determine the accurate 
properties of the oulcoming product and the optimum value to be extracted for a production run based on the available 
data. 

[0011] In one embodiment, including by way of example diesel fuel stocks with cetane improver or pour point im- 
so provers, a method is disclosed utilizing linear and non linear equations to more accurately predict the cetane number, 
pour point, and/or other properties of the resulting fuel product. 

[001 2] In another embodiment, tho present invention discloses a method of accurately predicting a characteristic of 
a product to be prepared by a processing plant, by providing an incoming material having at least one measured 
property, selecting at least one product having a desired characteristic to be prepared from said incoming materia! by 
55 said processing plant, selecting a process to be used by the processing plant to produce the product with its charac- 
teristic, and calculating a predicted value of the characteristic of the product utilizing a predictive equation. The pre- 
dictive equation is a validated non-linear equation generated by regressive analysis of an accumulation of data relating 
the measured property of the incoming material, the process, and the characteristic of the product. 
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[0013] As used hereinbelow, "reactor 0 should be construed in its broadest sense, to include mixing vessels, distillation 
columns, thermal cracking devices, etc. which may admix, treat, react or otherwise affect a material therein. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 

[0014] Figure 1 is a step block diagram of an exemplary production process. 
DETAILED DESCRIPTION 
'0 Pctrclcum Generally 

[0015] Crude oil. as a natural product, is very different in chemical composition depending upon the point of origin. 
Petroleum occurs throughout the world, and commercial fields have been located on every continent, crude oil from 
each region having very different properties. Petroleum is an extremely complex mixture and consists primarily of 
hydrocarbons as well as compounds containing nitrogen, oxygen, and sulfur. Most petroleums also contain minor 
amount of nickel and vanadium. Petroleums have a variety of different components with boiling ranges from about 20° 
C to above 650 a C. 

[0016] For an excellent discussion of the many and varied properties of crude oils, one may turn to, e.g. Marks' 
Standard Handbook for Mechanical Engineers, Tenth Edition,, T. Baumeister. McGraw Hill, 1 996 (ISBN 0-07-004997-1 ), 
"Petroleum and Other Liquid Fuels", J.G. Speight, pp. 7-10 through 7-14, which are incorporated herein by reference. 
[001 7] The chemical and physical properties of petroleum vary considerably because of the variations in composition. 
Crude oils are rarely used as fuel because the properties do not meet the specifications required by either furnaces, 
boilers, or engines. Crude oils, being complex mixtures of various complex chemicals, are not predictably mixed as 
the properties do not blend linearly. 

25 

Refining Processes 

[0018] Petroleum raw materials, e.g. crude oil from any of various sources in nature, are distilled into various fractions 
having differing commercial value. 

30 [001 9] Distillation separates the crude oil into fractions equivalent in boiling range to gasoline, kerosene, diesel fuel, 
lubricating oil, and residual. Thermal or catalytic cracking is used to convert kerosene, gas oil, or residual to gasoline, 
lower boiling fractions, and residual coke. Catalytic reforming, isomerization.alkylation, polymerization, hydrogenation. 
and combinations of these and other catalytic processes are used to upgrade the various refinery intermediates into 
improved gasoline stocks or distillates. The major finished products are usually blends of a number of stocks, plus 

35 additives. 

[0020] Typical end products include gasoline, jet fuel, diesel fuel, residual products, specialty products, and petroleum 
coke. 

[0021] Gasoline is a complex mixture of hydrocarbons that distills within the range of 100 to 400° F. Commercial 
gasolines are blends of straight-run, cracked, reformed, and natural gasolines. 
40 [0022] The specifications for gasoline (ASTM D439 and D48 1 4, which are incorporated herein by reference) provide 
for various volatility classes, varying from low-volatility gasolines to minimize vapor lock to high-volatility gasoline that 
permits easier starting during cold weather. 

[0023] Aviation gasoline has a narrower boiling range than motor gasoline. It has fewer low boiling and high boiling 
components. Aviation gasoline is defined by ASTM D910. 

[0024] Kerosene Is less volatile than gasoline and has a higher flash point, to provide greater safety in handling. 
Other quality tests are specific gravity, color, odor, distillation range, sulfur content, and burning quality. Specific tests 
for quality include flash point (minimum 115°F), distillation endpoint (maximum 572°F) sulfur (maximum 0.13 percent) 
and color (minimum +16) according to ASTM D187. 

[0025] Jet or Aviation turbine fuels are not limited by antiknock requirements and have wider boiling point ranges to 
assure greater availability for general aviation. Their properties are specified by ASTM D1655. Military requirements 
are as set forth in Military Specifications Mil-T-5624 (JP-4 and JP-5), Mil-T-83133 (JP-8), and Mil-P-87107 (JP-10). 
[0026] Diesel fuel is a liquid product distilled over the range of 350° C to 650°F. The carbon number ranges from 
about C10 to C18. The chemical composition of a typical diesel fuel and how it applies to the individual specifications 
• API gravity, distillation range, pour point, and flash point • are directly attributable to both the carbon number and the 
55 compound classes present in the finished fuel. 

[0027] Diesel fuels are measured in several aspects according to ASTM standards. They include API gravity (ASTM 
D1298), total sulfur (ASTM D2622), boiling point (ASTM D86), flash point (ASTM D93) s pour point (ASTM D97), hy- 
drogen content (ASTM D3701 ), cetane number (ASTM D6 1 3), acid number (ASTM D974), water and sediment content 
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(ASTM D1796), kinematic viscosity (ASTM D445), carbon residue on 10% residuum (ASTM D524), Ash (wt%) max 
(ASTM D482) and distillation temp (ASTM D86) 

[0028J The various diesel fuels for motor use require variability in performance since the engines range in size from 
small, high speed engines used in trucks and buses, to large, low-speed stationary engines for power plants. Thus, 
ASTM D975 provides for the specifications of a variety of diesel fuels. 

[0029] The combustion characteristics of diesel fuels are expressed in terms of the cetane number, a measure of 
ignition delay. A short delay (the time between injection and ignition) is desirable for a smooth running engine. Some 
diesel fuels contain cetane improvers, which usually are alkyl nitrates. The cetane number is determined by engine 
test (ASTM D61 3) or an approximate value, termed the cetane index, can be calculated for fuels which do not contain 
a cetane improver. 

[0030] The value of each of these components fluctuates on a daily basis, depending on supply and demand, market 
factors such as political disturbance in oil-producing regions, and weather, among others A lively futures market exists 
for various grades of crude oil and refined products. Keeping track of the cost of raw materials and prices of intermediate 
or finished products on a day-to-day basis is a daunting task. 

[0031] Processing crude oil is an extraordinarily complicated matter. The oil refinery separates the crude oil into 
individual compounds, or, more often, distillation fractions that consist of compounds with similar properties. 

Business Models and Solutions 

[0032] Operating a refinery or fuel blending and distillation plant is extraordinarily complex. Operating one in the 
most profitable manner possible is likewise exponentially more difficult. Several treatises are available to help in the 
understanding of the petroleum industry and refining, such as Petroleum Refinery Process Economics, R.E. Maples, 
1993, PennWell Publishing Company, Tulsa Oklahoma (ISBN 0-8781 4-384-X) and Petroleum Refining for the Non- 
technical Person, W.L. Leffler, 1979, 1985, PennWell Publishing Company, Tulsa Oklahoma (ISBN 0-87814-280-0), 
each of which is incorporated by reference. 

[0033] Multiple software programs are available to assist in the economic running of a refinery. They include several 
commercially available from Aspen Tech, Cambridge, Massachusetts 021 41 -2200. Aspen PIMS™ is a PC-based linear 
programming software module used by the petroleum and petrochemical industries. The software is capable of handling 
detailed operations planning, economic evaluation and scheduling activities based upon the cost and availability of 
raw materials, capacity considerations, and the demand for output. 

[0034] Aspen PIMS™ are a series of software tools for economic planning in the process industries. The system is 
designed to run on a Pentium™ class processor, or higher. Operating system platforms include Windows 95™, and 
Windows NT™. 

[0035] Aspen PIMS™ employs a linear programming (LP) technique utilizing a CPLEX™ optimizer, available from 
ILOG CPLEX Division. Incline Village, NV 89451, to optimize the operation and design of refineries, petrochemical and 
chemical plants or other facilities. It is stated to be useful for such varies processes as evaluation of alternative feed- 
stocks and product slates, optimization of operating decisions and product blending, and sizing of plant units in grass- 
roots and expansion studies. 

[0036] Aspen PIMS™ Scheduling Software assists in the preparation of detailed operating plans for material receipts, 
process operations, product blending and product shipments. 

[0037] SDPIMS™, also available from Aspen Tech, is software which models complex multiple-source, multiple- 
product, multiple-mode, multiple-tiered pricing, multiple-destination supply and distribution sytems and develops an 
optimized, least-cost solution for the entire network. Among other items, the time value of money is taken into account, 
as are time period lags for product shipments, and minimum, maximum, and target inventories. 
[0038] Aspen PIMS™ Refinery Scheduling System is software which is said to schedule refinery and petrochemical 
plant operations from feedstock arrivals to blendstock production, while Aspen PIMS™ Product Scheduling System is 
designed to schedule product blending and shipping activities . Aspen PIMS™ Pipeline Scheduling System is designed 
to schedule product deliveries through a pipeline network from product source through pumpstations and depots to 
final delivery point. 

[0039] In terms of the underlying crude oil related products and their properties, PassMan™ also available from 
Aspen, is a PIMS™ crude oil assay manager whose function is to serve as a manager of a crude oil library and to 
output a table of crude oil data that is recut to the needs of the end user either for input into their linear production 
model or other applications. Other assay databases are available. 

Optimization 

[0040] Sunset Software Technology™, based in San Marino, California 91108, supplies multiple optimization algo- 
rithm based products. Sunset's products include linear, binary, mixed-integer, interior point (barrier) and quadratic 
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programming products and services, which operate on platforms ranging from PC's to UNIX workstations. They include 
products currently marketed under the following names: XA Linear Optimizer System™, XA Binary and Mixed Integer 
Solver™, XA Barrier Solver™, XA Quadratic Solver™, XA Parallel MIP Solver™, and XA Callable Library™. 
[0041] CPLEX™ Base Development system, available from ILOG, includes several commercial optimizers - primal 

5 Simplex, dual Simplex, and network Simplex solvers for linear programming problems. This program is in one embod- 
iment configured with the CPLEX Callable Library™ , which provides CPLEX algorithms in a library of CPLEX algorithmic 
and utility routines. The CPLEX Mixed Integer Solver Option™ includes the capability to solve problems with mixed 
integer variables (general or binary), utilizing algorithms and techniques, including cuts (cliques & covers), heuristics, 
and a variety of branching and node selection strategies. CPLEX Barrier/QP Solver™ is a primal-dual log barrier al- 

'O gorithm with predictor corrector said to be useful for solving certain classes of linear programming models and quadratic 
programming problems. 

[0042] AIMMS™ (Advanced Integrated Multidimensional Modeling Software), available from Paragon Decision Tech- 
nology B.V , P.O. Box 3277, 2001 DG Haarlem, The Netherlands, is a software package which allows modelers to 
create functional analytic decision support applications. AIMMS contains a graphical model explorer, which builds and 
*s maintains complex (optimization-based) modeling applications. AIMMS is said to be capable of modeling a particular 
(optimization-based) decision support problem, creating an end-user interface around the model suitable for use by 
end users. 

[0043] AMPL™ modeling language, available from CPLEX (a division of ILOG) Is an algebraic modeling language 
for linear, nonlinear, and integer programming problems. It is said to be useful for optimization model types including 
20 linear programming problems, network problems, mixed integer programming problems, quadratic programming prob- 
lems, and general non-linear programming problems. 

[0044] The General Algebraic Modeling System (GAMS™), available from GAMS Development Corporation, 1217 
Potomac Street NW, Washington, DC 20007 USA is software capable of modeling linear, nonlinear and mixed integer 
optimization problems. GAMS can solve LP, MIP and different forms of NLP models. 
25 [0045] MPL™ (Mathematical Programming Language), available from Maximal Software, Inc., 2111 Wilson Boule- 
vard, Suite 700, Arlington, VA 22201, U.S.A. is a modeling system that permits the construction of complex models, 
involving thousands of constraints which allows the import of data directly from a database and then export of the 
solution back into the database. 

[0046] Models developed in MPL can be used with nearly all LP-solvers on the market today as MPL supports a 

30 number of industrial strength solvers. 

[0047] The mathematical technique known as linear programming is commonly used by many of the above programs 
to solve a variety of industrial and scientific problems by arriving at an "optimal solution". Linear programming ("LP") 
has existed from about the 1940's. It works by creating an LP "model" which represents some situation that is then 
solved to discover the optimum plan. 

35 [0048] A valid LP model must have four elements. First, there must be an objective function. Generally for business 
models the value to maximize is profit, and the value to be minimized is usually cost or distance. Each activity in the 
model contributes to this objective, either favorably or unfavorably. Second, there are limited resources. For example, 
a machine can only run for 24 hours in a day, or only so much material is available to purchase. Third, there must be 
linear relationships between activities and these resources. For example, one relationship might be the number of 

40 machine hours which should be operated and how much material should be bought. Finally, there must be an assump- 
tion of certainty, i.e. an assumption that these conditions in the model will be resolved feasibly. A more sophisticated 
technique known as stochastic programming exists to handle probability-based programming. 
[0049] A variety of industrial LP applications have been developed to solve varying requirements. For example, in 
the field of product planning, one can plan an appropriate mix by solving the LP for optimal production quantities of 

«5 products subject to resource capacity and market demand constraints. For blends, one can solve for optimal proportions 
of ingredients for products such as gasoline, foods, livestock feeds, subject to certain minimal requirements. 
[0050] For distribution, one can use LP to solve for optimal shipping assignments from factories to distribution centers 
or from warehouses to retailers. For location planning, e.g. of facilities, one can determine the optimal location of a 
plant or warehouse with respect to total transportation costs between various alternative locations and existing supply 

so and demand sources. 

[0051] For process control, one can use LP models to, e.g.solve for the cutting pattern that minimizes the amount 
of scrap material, given the dimensions of a roll or sheet of raw material. For scheduling, one can use LP to determine 
the minimum-cost assignment of workers to shifts subject to varying demand. For vehicles, one can assign available 
vehicles to jobs and determine the number of trips to make, subject to vehicle size, availability, and demand constraints. 
55 Similarly, for routing, one can solve for the optimal routing of a product through a number of sequential processes, 
each with its own capacities and characteristics. 

[0052] For production planning, one can solve for minimum-cost production scheduling for a given work force, taking 
into account inventory carrying and subcontracting costs. The management of production and work force may be 
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accomplished by LP by solving for minimum-cost production scheduling, taking into account hiring and layoff costs as 
well as inventory carrying, overtime, and subcontracting costs, subject to various capacity and policy constraints. Fur- 
thermore, one can solve for optimal staffing for various categories of workers, subject to various demand and policy 
constraints 

[0053] More sophisticated techniques involve the modification of standard LP techniques, but which relax some of 
the assumptions of the basic LP model. 

[0054] Integer, binary, and mixed integer program modeling allow for activities that may only be conducted incre- 
mentally. With shipping for example, if one truckload is shipped, costs of the truck must be calculated whether one 
pallet or 24 pallets are shipped. 

[0055] Non-linear programming allows for non-linear relationships between activities and constraints, while stochas- 
tic programming allows for uncertainty. 

[0056] A presently preferred optimizer program for the instant invention is What'sBest!™ 4.0, available from Lindo 
Systems. Inc.. 1415 North Dayton Street, Chicago, IL 60622, USA, which is an add-in to Excel™ (available from Mi- 
crosoft. Inc. Redwood. WA) that allows the building of large scale optimization models in a free form layout within a 
spreadsheet. What'sBest!' M combines the linear, nonlinear and integer optimization with Microsoft™ Excel 1M . 

Optimization in the Petroleum Industry 

[0057] As noted above., there are many uses for various petroleum products, and each end use, e.g of diesel fuel 
has vatying requirements for the product for safe, environmentally conscious, and economical use. 
[0058] However, the fuel stocks on hand may frequently be limited to nonoptimal fuel stocks. For example, the refiner 
may need to determine whether it is more economical to blend a distillate component, such as Light Cycle Oil, into 
diesel or into residual product. Each of the fuels has a different value, and each of the raw materials (diesel fuel stocks, 
residual stocks, cctanc improvers, etc) has a different value as a blending component. 

[0059] The traditional approach has included mixing stocks in a blending operation to attempt to solve supply issues, 
but the accurate prediction of non-linear properties combined with a maximization of profit potential has been limited 
by the linear nature of the optimization software. There has been until now no validated way to accurately predict the 
properties of complex systems, e.g. blended fuels, additive enhanced fuels, and the like. In production, refiners are 
frequently forced to remeasure the data and compare it to a predicted value, making production decisions based upon 
these predicted properties delayed and/or flawed. 

[0060] In one embodiment of the present invention, a complex production process can be more effectively managed 
from a business standpoint by assessing the relative values of incoming materials or components of the final product, 
and accurately predicting the properties of the final products. 

[0061 ] Turning to Figure 1 , a hypothetical production scenario is represented by the step block diagram. Various raw 
material sources, each having different properties, are labeled A, B, C, D and E. Production reactors R1 and R2 are 
available for the mixing or other processing of the raw materials. Oesired end or final products are I, II and III. 
[0062] Raw materials A-E may be virtually any raw material - crude oils for distillation, gasoline stocks for blending, 
additives and fuels for specialty fuel production, sands (e.g. silica) for raw glass or ceramic production, paint pigments 
and solvents for paint systems, particulates for alloy manufacture by tape casting, die casting, sintering, annealing, 
grains for cereal or bread production, nutraceuticals for vitamin manufacture, etc. The raw materials may be combinable 
by mixing, reacting, or otherwise commingling; or may be separatable by fractionation, distilling, cracking, or the like. 
[0063] Each applicable industry to which this invention may be applied will have its own known methods of mixing, 
producing, refining, etc.. which result in known products. The following examples, while most exemplary of the invention 
in the petroleum blending and additive industry, should not be construed as limiting, as the invention has far-ranging 
uses which one skilled in the art, having regard for this disclosure, will easily be able to achieve. 
[0064] Fuels are blended differently for seasons and uses, as illustrated above. There are several differentcategories 
of blendstocks available for fuel use - aviation, kerosene, gasoline, diesel, and residual are the most common. 
[0065] The blending of fuels is not represented by a linear relationship. Given the complex nature of the various 
petrochemicals present in fuelstocks, simply mixing a fuel component of a known cetane number, for example, with 
another of known cetane number does not necessarily result in an intermediate cetane number fuel. It is possible for 
a blended fuel to have a lower cetane number than either of the component diesel fuels, resulting in a loss in value 
(negative blending). Thus, it is vital to the economic survival of the refiner or manufacturer to have an accurate prediction 
of the properties and values of the resulting product. 

[0066] By way of nonlimiting example, one possible manufacturing operation is the creation of a diesel fuel, by blend- 
ing fuel streams and/or by the inclusion of a cetane improver. The standard method to measure the cetane number is 
the use of a cetane engine as described in ASTM D613. It is possible presently to predict cetane number with the 
cetane index equations indicated in ASTM D976 and ASTM D4737. However, such cetane index prediction models 
generally have an error of at least 2 cetane numbers. This is a costly error which can result in "cetane giveaway". Not 
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only is cetane index a poor predictor, it does not take into account the improvement in cetane number caused by 
additives. 

[0067] Thus, if a pipeline has a requirement of a cetane number of 40. the cetane index must be at least 42 to assure 
the minimal requirements of the pipeline are being met. Table I, below indicates the cetane standards for the Colonial 
Pipeline Company, of Atlanta, Georgia. 



TABLE I 



Colonial Pipeline Company 

Specifications for Fungible Low Sulfur Diesel Fuel (cetane only) 
Issue No. 3; Grade 74 




ASTM Test 


Test Results 


PRODUCT 
PROPERTY 


Method 


Min 


Max 


Notes 




Notes 












Cetane Number 


D613 




40 






Cetane Index 


D976 




42 




4 


NOTES 













(4) Where celane number by test method D613 is not available, test method D976 or D4737 can be used as an approximation. Minimum cetane 
index of 42 accounts lor +/- 2 accuracy of the cetane index methods in approximating cetane number. 



[0068] As may be seen from the above, one supplying fuel to the Colonial Pipeline using the standard Cetane Index 
calculation must supply fuel of a minimum of 42 cetane index, in effect giving away the commercial value of 38, and 
potentially 42 cetane number diesel fuel. Cetane Number must be determined using a Cetane Engine after the blend 
30 has been formulated and prepared, which is often too late to be of substantial economic advantage during preparation 
and blending. 

[0069] However, with the novel invention, sufficient data is collected on a wide range of properties relating to a 
plurality of fuel components such that an extremely accurate prediction of the cetane number may be made. 
[0070] Accurate prediction of the properties of the resultant product is important to assist in optimizing the output of 
35 a plant which produces such products; for example, the operator of such a plant may make more of a higher priced 
product than a lower priced product and thereby maximize profit, so long as the operator is assured that the final 
product will meet specifications. By means of this invention, additional components or reblending after formulation is 
not needed. 

[0071] ladomrnatos and Goacher in "Equations for predicting the cetane number of diesel fuels from their physical 
40 properties", Fuel, Vol. 74 No. 7, pp. 1 083-1 093 (1 995) derived twenty-two equations for predicting the cetane number 
of diesel fuels. Likewise, Maxwell et all in "How to accurately predict cetane numbers of diesel-fuel blend stocks",0// 
and Gas Journal, November 3, 1 969 developed predictive equations for cetane value. These articles are incorporated 
by reference as if fully set forth herein. 

[0072] Each of these articles, however, fails to provide a reliable prediction of the variety of properties required by a 
45 complex system of physical parameters. They rely primarily upon the cetane index, which alone is notoriously inaccu- 
rate. 

[0073] An exemplary prediction of the cetane number of diesel fuel containing additives will be discussed as follows. 
[0074] Turning now to Table II, several diesel fuel properties are indicated which are important for meeting the various 
ASTM or military specifications for such fuels (see, e.g. ASTM D975). 

50 

TABLE II- 

DIESEL FUEL PHYSICAL PROPERTIES 

Hydrocarbon Type 
Carbon content 

55 Hydrogen content 

Sulfur content 
Nitrogen content 
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TABLE II - (continued) 
DIESEL FUEL PHYSICAL PROPERTIES 

API (specific) gravity 
5 Distillation range 

Cetane number 
Aniline point 
Heat content 
Kinematic viscosity 
10 Cloud point 

Pour point 
Flash point 



[0075] Data were collected on 1 54 low sulfur diesel fuels to derive cetane-prediction equations for diesel fuels with 
and without the addition of a cetane improver The equations were derived by the use of statistical analysis including 
multiple linear regression to derive equations to predict the cetane number of fuels not in the original data set Equations 
containing different combinations of variables were developed to accommodate laboratories with different analytical 
capabilities. The invention thus provides prediction models using from 4 to 8 or more input variables. 
[0076J The equations were then validated with twenty new fuels not in the original data set of 154 fuels. 
[0077] The equations derived are a set of non-linear calculations as follows: 
Definitions: 

CN = cetane number (ASTM D613) 

CN 0, 1000, 2500, 5000, 7500 ! 10000 = cetane number with ppmv of cetane improver DII-3™ 

T90= temperature at which 90% of the material boiled off (°C) 

T50= temperature at which 50% of the material boiled off (°C) 

T10= temperature at which 10% of the material boiled off (°C) 

ANPT = aniline point (°C) 

VISC40 = viscosity at 40° C (ASTM D445) 

D976 = cetane index as provided by ASTM D976 

CLOUD = temperature at which wax crystallization occurs {'C) 

SPGR = specific gravity (no dimension) 

FIAAROM = aromatic content (vol. %) measured by ASTM D1319 

[0078] The exemplary preferred models for diesel fuel cetane number prediction upon inclusion of a cetane improver 
follow: 

MODEL I (wherein there are 5 inputs) 
[0079] 

CN 0 = T90 x (-0.073074) + ANPT x 0.307979 + VISC40 x (-1 .1 521 77) + D976 x 
0.49566 + CLOUD x 0.183557 + 31 .572661 



CN 1000 = T90 x (-0.093739) + ANPT x 0.273761 + VISC40 x (-0.951212) + D976 x 
50 0.588985 + CLOUD x 0.206887 + 40.162203 

CN 2500 = T90 x (-0.096809) 4- ANPT x 0.2521 26 + VISC40 x (-0.988991 ) + D976 x 
55 0.670028 + CLOUD x 0.194338 + 42.384571 
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CN 5000 = T90 X (-0.090088) + ANPT x 0.308413 + VISC40 x (-1.183015) + D976 x 
0.6461 13 + CLOUD x 0 16953 + 41 .297981 

CN 7500 = T90 X (-0.07601 5) + ANPT x 0.362383 + VISC40 x (-1 .306567) + D976 x 
0.649497 + CLOUD x 0.137109 + 35.771915 

CN 10000 « T90 x (-0.069705) + ANPT x 0.41 8994 + VISC40 x (-1 .257463) + D976 x 
0.603263 + CLOUD x 0.1 08566 + 33.746498 



MODEL 2(5 inputs) 
[0080] 

20 

CNO = T50 x 0.170022 + T90 x (-0.090224) + SPGR x (-238.35652) + FIIAROM x 
(-0.119872) + CLOUD x 0.218696 + 237776061 

CN1000 = T50 x 0.180185 + T90 x (-0.106238) + SPGR x (-256.07975) + FIIAROM x 
(-0.112265) + CLOUD x 0.241931 + 259.828051 

CN25O0 = T50 X 0.186051 + T90 X (-0.107282) + SPGR x (-264.7448) + FIIAROM X 
(-0.121278) + CLOUD x 0.230962 + 269.91 1717 

CN5000 = T50 x 0.192B09 + T90 x (-0.103698) + SPGR X (-267.85828) + FIIAROM x 
(-0.148856) + CLOUD x 0.21 1281 + 273.437865 

CN7500 = T50 x 0.207947 + T90 x (-0.092333) + SPGR x (-283.78141) + FIIAROM x 
(-0.169967) + CLOUD x 0.18288 + 282.147518 

CN10000 = T50 x 0.216317 +T90 x (-0.086368) + SPGR x (-285.49607) + FIIAROM x 
(-0.193056) + CLOUD x 0.157395 + 281.674744 

so MODEL 3 (5 inputs) 
[0081] 

55 CNO = T10 X 0.052862 + T50 x 0.132853 + T90 x (-0.073079) + SPGR x (-319.241664) 

+ CLOUD X 0.1 98433 + 295.090575 
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CN1000 = T10 x 0.055293 + T50 x 0.1 39078 + T90 X (-0.08731 8) + SPGR x 
(-332.515775) + CLOUD x 0.221971 + 313.561752 

CN2500 = T10 x 0.065233 + T50 x 0.135657 + T90 x (-0.0841 19) + SPGR x 
(-347.966881) + CLOUD x 0.208466 + 328.012419 

CN5000 = T10 x 0.091 782 + T50 x 0.1 18204 + T90 x (-0.06947) + SPGR x 
(-371.389545) + CLOUD x 0.18161 +344.864624 

CN7500 = T10 x 0.120328 + T50 x 0.105859 + T90 x (-0.045566) + SPGR x 
(-403.830741) + CLOUD x 0.146444 + 363.855067 

CN10000 = T10 x 0.137699 + T50 x 0.099245 + T90 x (-0.03274) + SPGR x 
(-421.974629) + CLOUD x 0.115835 + 374.491796 



MODEL 4 (3 Inputs) 
[0082] 



CN0 = ANPT X 0.281 1 22 + VISC40 X (-1 .0301 39) + D976 X 0.651 89 
CN1000 = ANPTx 0.240331 + VISC40 x (-0.909077) + D976 x 0.794587 
CN2500 = ANPT x 0.217108 + VISC40 x (-1 .074988) + D976 x 0.914676 
CN5000 = ANPT x 0.2631 56 + VISC40 x (-1 .3351 08) + D976 x 0.936531 
CN7500 = ANPT x 0.30481 8 + VISC40 X (-1 .322044) + D976 X 0.934638 
CN10000 = ANPT x 0.3534 + VISC40 x (-1 .25902) + D976 x 0.905473 



MODEL 5 (7 inputs) 
so [0083] 

CN0 = T1 0 x -0.565869 + T50 x T90 x 0.000458 + (T90) 2 x (-0.000992) + T1 0/SPGR 
55 X 0.526071 -f T90/SPGR x 0.426314 - 64.399604 

CN1000 = T10 x -0.443266 + T50 x T90 x 0.000483 + (T90) 2 x (-0.001215) + T10/SPGR 



10 



BNSDOCfD: <EP 1 1C21B7A2_I_> 



EP 1 102 187 A2 

x 0.425422 + T90/SPGR x 0.536406 - 81 .532783 



CN2500 = T10 x -0.473346 + T50 x T90 X 0.000468 + (T90) 2 x (-O.001226) + T10/SPGR 
x 0.458768 + T90/SPGR x 0.542029 - 79.539029 



CN5000 = T10 x -0.582345 - T50 X T90 X 0.000399 + (T90) 2 X (-0.001 13) + T10/SPGR 
x 0.573578 +T90/SPGR x 0.50863 - 73.521591 



CN7500 = T1 0 x -0.62471 3 + T50 x T90 x 0.000359 + (T90) 2 x (-0.001 1 25) + T10/SPGR 
x 0.632041 + T90/SPGR x 0.524471 - 79.914841 



20 CN10000 = T10 x -0.722336 4- T50 x T90 x 0.000346 + (T90) 2 x (-0.001046) + 

T10/SPGR x 0.727228 + T90/SPGR x 0.487577 - 74.343479 

[0084] Cetane numbers at intermediate concentrations of cetane improver are calculated by interpolating the values 
25 from these formulae. The formulae may be reduced to lower numbers of significant figures, e.g. 2, 3, or 4 with minimal 
loss of accuracy of prediction 

[0085] This method allows the calculation of cetane numbers of the product from the properties of the components, 
and does not rely on the cetanes. As discussed infra, many of the ingredients in a complex chemical system react in 
a nonlinear fashion and result in unexpected properties. Accumulation of data and the subsequent inclusion of correc- 
30 tions in a nonlinear model enables the inventive method to provide a far more accurate prediction of the properties of 
the resultant composition or system. 

[0086] The formulas which follow, containing fewer significant digits, also provide the benefits of the instant invention 
Model 1 

35 

[0087] 

CN0 = -0.07307 x T10 + 0.3080 x ANPT- 1 .152 x VISC40 + 0.4957 x D976 + 
0.1836 x CLOUD +31.57 



CN1000 = -0.09374 x T10 + 0.2738 x ANPT - 0.9512 x VISC40 + 0.5890 x D976 + 
0.2069 x CLOUD + 40.16 



CN2500 = -0.09681 x T10 + 0.2521 x ANPT - 0.9890 x VISC40 + 0.6700 x D976 + 
0.1943 x CLOUD + 42.38 



CN5000 = -0.09009 x T10 + 0 3084 x ANPT - 1 .183 x VISC40 + 0.6461 x D976 + 
0.1695 x CLOUD + 41.30 
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CN7500 = -0.07602 x T10 + 0.3624 x ANPT - 1 .307 x VISC40 + 0.6495 x D976 + 
0.1371 x CLOUD + 35.77 

5 

CN10000 = -0.06971 x T10 + 0.4190 x ANPT - 1 .257 x VISC40 + 0.6033 x D976 + 
0.1086 x CLOUD + 33.75 

W 

[0088] A generic formula is also provided which is suitable for the practice of the instant invention: 
CN(X) = T10 x (-0.06971 to -0.09681 ) + ANPT x (0.2521 to 0.4190) + VISC40 x 

15 

(-0.9512 to -1 .307) + D976 x (0.4957 to 0.6700) + CLOUD x (0.1086 to 0.2069) + (31.57 to 42.38), 

where (X) equals ppmv of cetane improver. 
20 Model 2 
[0089] 
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CNO = 0.1 700 x T50 - 0.09022 x T90 - 238.4 x SPGR - 0.1 1 99 x FIAAROM + 
0.21 87 x CLOUD +237.8 

CN1000 = 0.1802 X T50 - 0.1062 X T90 - 256.1 x SPGR - 0.1 123 X FIAAROM + 
0.2419 x CLOUD + 259.8 

CN2500 = 0.1861 x T50 - 0.1073 x T90- 264.7 x SPGR - 0.1213 x FIAAROM + 
0.2310 x CLOUD + 269.9 

CN5000 = 0.1928 x T50 - 0.1037 x T90 - 267.9 x SPGR - 0.1489 x FIAAROM + 
0.2113 x CLOUD + 273.4 

CN7500 = 0.2079 x T50 - 0.09233 x T90 - 283.8 x SPGR - 0.1 700 x FIAAROM + 
0.1829 x CLOUD + 282.1 

CN10000 = 0.2163 x T50 - 0.08637 x T90 - 285.5 x SPGR - 0.1931 x FIAAROM + 
0.1574 x CLOUD + 281.7 



[0090] A generic formula is also provided which is suitable for the practice of the instant invention: 

55 

CN(X) = T50 x (0.1700 to 0.2163) + T90 x (-0 08637 to -0.1073) + SPGR X (-238.4 
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to -285.5) + FIAAROM x (-0.1 123 to -0.1931) + CLOUD x (0.1574 to 0.2419) +(237.8 to 282.1), 
where (X) equals ppmv of cetane improver. 

5 

Mode! 3 
[0091] 

w 

CN0 = 0.05286 x T10 + 0.1329 x T50 - 0.07308 x T90 - 319 2 x SPGR + 0.1984 x 

CLOUD + 295.1 

CN1000 = 0.05529 x T10 + 0.1391 x T50 - 0.08732 x T90 - 332.5 x SPGR + 
0.2220 x CLOUD + 313.6 

CN2500 = 0.06523 x T1 0 + 0. 1 357 x T50 - 0.0841 2 x T90 - 348.0 x SPGR + 
0.2085 x CLOUD + 328.0 

CN5000 = 0.091 78 x T1 0 + 0. 1 1 82 x T50 - 0.06947 x T90 - 371 .4 x SPGR + 
0.1817 x CLOUD + 344.9 

CN7500 = 0.1203 X T10 + 0.1059 x T50 - 0.04557 X T90 - 403.8 X SPGR + 0.1464 x 

CLOUD + 363.9 

CN1000 = 0.1 377 x T1 0 + 0.09925 x T50 - 0.03274 x T90 - 422.0 x SPGR + 
0.1 158 x CLOUD + 374.5 

40 [0092] A generic formula is also provided which is suitable lor the practice ol the instant invention: 

CN(X) = T10 x (0.05286 to 0.1377) + T50 x (0.09925 to 0.1391) + T90 (-0.03274 to 
45 -0.08732) + SPGR x (-319.2 to -422.0) + CLOUD x (0.1 158 to 0.2220) + (295.1 to 374.5), 

where (X) equals ppmv of cetane improver. 
Model 4 

so 

[0093] 

CN0 = 0.2811 xANPT- 1 .030 x VISC40 + 0.6519 x D976 

55 

CN1 COO = 0.2403 x ANPT - 0.9091 X VISC40 + 0.7946 X D976 
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CN2500 = 0.21 71 x ANPT - 1 .075 x VISC40 + 0.91 47 x 0976 

5 CN5000 = 0.2632 x ANPT • 1 .335 x VISC40 + 0.9365 X D976 

CN7500 = 0.3048 x ANPT - 1 322 x VISC40 + 0.9346 x D976 

10 

CN10000 = 0.3534 x ANPT - 1 .259 x VISC40 + 0.9055 x D976 
[0094] A generic formula is also provided which is suitable for the practice of the instant invention. 

15 

CN(X) = ANPT x (0.2171 to 0.3534) + VISC40 (-0.9091 to -1.335) + D976 x (0.6519 to 0.9365), 

where (X) equals ppmv of cetane improver. 
20 Model 5 
[0095] 

25 CN0 = -0.5659 x T1 0 4- 0.000458 x T50 x T90 - 0.000992 x (T90) 2 + 0.5261 x 

T10/SPGR + 0.4263 x T90/SPGR- 64.40 

30 CN1 000 = -0.4433 x T10 + 0.000483 x T50 x T90 - 0.001215 X (T90) 2 + 0.4254 x 

T1 0/SPGR + 0.5364 X T90/SPGR - 81 .53 

35 CN2500 = -0.4733 X T1 0 + 0.000468 X T50 X T90 - 0.001 226 X (T90) 2 + 0.4588 X 

T1 0/SPGR + 0.5420 x T9 0/SPGR - 79.54 

40 CN5000 = -0.5823 x T10 + 0.000399 x T50 x T90 - 0.001 130 x (T90) 2 + 0.5736 x 

T1 0/SPGR + 0.5086 x T90/SPGR- 73.52 



CN7500 = -0.6247 x T10 + 0.000359 x T50 x T90 - 0.001 125 x (T90) 2 + 0.6320 x 
T1 0/SPGR + 0.5245 x T90/SPGR- 79.91 



50 CN10000 = -0.7223 x T10 + 0.000346 x T50 x T90 - 0.001046 x (T90) 2 + 0.7272 x 

T1 0/SPGR 4- 0.4876 x T90/SPGR - 74.34 

[0096] A generic formula is also provided which is suitable for the practice of the instant invention: 

55 

CN(X) = T10 x (-0.4433 to -0.7223) 4- T50 x T90 x (0.000346 to 0.000483) + (T90) 2 x 
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(-0.000992 to -0.001226) + T10/SPGR x (0.4254 to 0.7272) -f T90/SPGR x (0.4263 to 0.5420) + (-64.40 to -81.53), 
where (X) equals ppmv of ceiane improver. 

[0097] Turning now to Table III below, a comparison of the values of the cetane number of 20 fuels, calculated by 
the prior art methods (D-976, 0-4737, and CGSB (Canadian General Standards Bureau method) versus the instant 
models is illustrated. 



TABLE ill 



w 



20 



30 



Fuel 


CN 


D-976 


D-4737 


CGSB 


Model 1 


Model 2 


Model 3 


Model 4 


Moripl ^ 

IVIVUCI >J 


A 

A 


AO C 

4B.D 


49.8 


51.6 


50.1 


48.1 


46.9 


48.4 


48.3 


48.7 


r> 
D 


47. 8 


48.2 


48.0 


46.5 


45.7 


45.1 


45.2 


46.4 


46.0 


C 


41 .2 


43.7 


43.1 


40.7 


41 .3 


41 .4 


41.3 


41.1 


41.2 


D 


45.0 


47.0 


47.4 


41 .0 


43.8 


46.8 


46.1 


43.7 


45.9 


t 




46.1 


45.7 


44.5 


44.0 


43.5 


43.4 


44.7 


44.6 


F 


44.9 


45.2 


45.2 


44.4 


45,2 


45.4 


45.0 


44.4 


44.5 


G 


42.1 


48.2 


48.0 


46.4 


46.2 


46.0 


45.6 


46.6 


46.0 


H 


47.6 


51.7 


52.9 


52.6 


51.6 


50.9 


51.0 


51.2 


50.8 


1 


40.2 


42.0 


40.7 


38.7 


40.2 


40.2 


39.7 


39.3 


39.0 


J 


43.2 


43.4 


42.4 


40.4 


41.9 


42.1 


41.5 


41.0 


40.7 


K 


51.7 


55.8 


55.7 


52.8 


53.1 


53.3 


52.3 


53.9 


52.6 


L 


43.6 


46.8 


47.0 


43.7 


44.7 


45.1 


45.2 


44.7 


45.4 


M 


37.0 


41.5 


40.9 


40.7 


37.7 


37.3 


36.8 


39.6 


39.3 


N 


46.0 


47.0 


47.1 


44.0 


44.5 


44.7 


44.9 


44.8 


45.4 


O 


47.9 


50.7 


51.4 


49.2 


48.1 


48.4 


48.2 


48.7 


48.7 


P 


43.4 


45.2 


45.4 


44.2 


45.2 


45.8 


45.0 


44.3 


44.5 


Q 


49.0 


49.6 


50.4 


49.6 


47.8 


47.8 


47.3 


48.4 


48.0 


R 


52.6 


52.2 


54.1 


53.2 


48.2 


48.1 


48.0 


50.9 


50.9 


S 


40.8 


42.2 


41.3 


42.0 j 


40.0 


39.1 


38.9 


41.0 


40.6 


T 


47.2 


48.5 


48.2 


45.8 


46.2 


45.6 


46.0 


46.6 


46.5 


RMSE 




2.6 


2.8 


2.3 


1 .9 


2.0 


1.9 


1.8 


1.6 



[0098] It is clearly seen that the root mean square error (RMSE) of the inventive method is significantly below that 
of any of the prior art methods for predicting the cetane number of diesel fuel with a cetane improver therein. 
[0099] In a preferred embodiment of the present invention, there are at least 8 component properties which are 

40 collected from actual samples, and resulting blended products. The properties are then entered into a database or 
spreadsheet and correlated to ceiane number by linear regression analysis. In another preferred embodiment, there 
are at least 4 (four) properties measured for each component used in the final blend. Five, six, seven, or more properties 
are also acceptable, and there may be as many as required to formulate a statistically proper prediction equation. 
[0100] Preferably, the prediction calculation is tied into a spreadsheet, such as Excel™, even more preferably along 

45 with additional calculations such as cost. 

Initially, the product lo be produced, and its characteristics are selected. Such characteristics include typically such 
qualities as minimum and maximum volumes, price of the end product, minimum and maximum API. maximum sulfur 
content, minimum flash point, maximum aromatics content, maximum freeze point, T90, cloud point, pour point, vis- 
cosity at 100 degrees C (ASTM D445), viscosity at 40 degrees C, vanadium content, and cetane number. 

so [0101] Then, the appropriate feedstock to create the desired product is provided to the model. Characteristics of the 
feedstock to be selected for fuel production include amount available, purchase price, API, sulfur content, aromatics, 
naphthalene content, smoke, vanadium content, aniline point, cloud point, flash point, freeze point, pour point, viscosity 
at 100°C, viscosity at 40°C, distillation T10, distillation T50, distillation T90, and cetane number. 
[0102] Then, the product requirements and raw materials data are inputted into the database either manually or by 

55 importation from another database in a known manner. Table IV illustrates an exemplary entry table for cetane require- 
ments for a finished product. 
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TABLE IV - 



10 



15 



20 



25 



30 



35 



Cetane Specification Input Table for Diesel Fuel 


CETANE 


CALCULATIONS 




Data Required: 


Method 1 


Ethyl 


API, T50, 


T90, 


An PI, Vise, 


Cloud 


Method 2 


Ethyl 


API, T50, 


T90, 




Cloud, Arom 


Method 3 


Ethyl 


API, T50, T10, 


T90, 




Cloud 


Method 4 


Ethyl 


API, T50, 




AnPt, Vise 




Method 5 


Ethyl 


API.T50.T10, 


T90 






ASTM Method 


D-976 


API, T50 









[0103] Depending on the data available, the appropriate method for calculating the required input ingredients is 
selected. Methods 1-5 are according to the present invention for calculating cetane requirements, the ASTM method 
D-976 is according to the prior art. It is notable that D-976 requires only two inputs. API gravity and T50. 
[01 04] If additional requirements are desired, they may be inputted via a similar spreadsheet table. Any other desired 
information, including processing information, costs, times, volumes, etc. may be included. An additionally preferred 
embodiment includes calculations for cold flow improver for diesel fuel once the variables have been calculated utilizing 
the multiple regression analysis of the instant invention. 

[0105] Once the product characteristics have been inputted, a computing device (preferably a Pentium based PC 
running a form of Windows or equivalent software) then may calculate the required feedstock to manufacture the 
product. The calculations may be run in reverse, i.e. to determine what may be produced from the feedstock on hand. 
[01 06] Furthermore, the instant invention allows for the "tuning" of the equations for particular refineries, fuel sources, 
geographical variances, etc. with the collection of additional data and verification through actual samples. Each of the 
variables may be adjusted by one of skill in the art having regard for this disclosure by regressive analysis or similar 
means. 

[0107] The efficient and profitable operation of a large scale manufacturing facility is illustrated as follows. A target 
cetane number for fuel is selected for a particular application, e.g. 45. 

[01 08] The fuel may be made a number of ways, by mixing various components, and/or by adding a cetane improver 
to the blend. Turning to Table V. which follows, the following pricing information is assumed to be known (not actual 
amounts): 



40 



Low Sulfur Diesel market price $23.10 / bbl 
Residual Fuel Oil market price $16.38/bbl 
D1I-3 market price $0.75/pound 

Nominal refinery diesel output: 50,000 barrels per day (B/D) 



[0109] The method according to the instant invention determines the precise amounts of components of the types 
of Virgin Distillate, Light Cycle Oil (LCO), and Cetane Improver to add to reach the desired property of 45 cetane 
number. However, profit optimization can help select the more economical production route, when ordinarily one such 
45 route might have been ignored. Case One is a traditional blend to result in the desired CM, and requires that a large 
amount of LCO be placed into the less profitable residual product. Whereas, in Case Two, the inclusion of Dll-3™ 
cetane improver allows more LCO to be incorporated into the Diesel Pool, and less LCO into the less profitable Resid 
Pool. The sums saved are significant. 
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Product data: 



Diesel price: 523.10/bbl 
Resid price: $16.38/bbl 
Target CN: 45 



Case One, low usage rate of inexpensive coraponet LCO 



Component 


Bbl placed in Diesel 
Diesel pool Rev, $ 


Bbl placed Resid 
in 

Resid pool Rev, $ 


Total 
Rev, S 


Virgin distillate 
LCO 


35,000 808,500 
5,000 115,500 


0 0 
10,000 163,800 


808,500 
279,300 


Total Rev: 1,087,800 



Case Two, high usage rate of inexpensive componet LCO 



Component 


Bbl placed in Diesel 
Diesel pool Rev, $ 


Bbl placed Resid 
in 

Resid pool Rev, S 


Total 
Rev, $ 


Virgin distillate 

LCO 

DII-3 


34,950 807,345 
9,950 229,845 
100 -25,295 


0 0 
5.000 81,900 


807,345 
311,745 
-25,295 


Total Rev: 1,093,796 



Results (comparing Case Two to Case One): 

Savings per bbl: $0.12 
Savings per day: $5,996 
Savings peryr: $2,188,358 



[0110] The entire disclosure ol all applications, patents, available software, ASTM and other standards, and publi- 
cations cited above and throughout this application are hereby incorporated by reference. 

[0111] The preceding examples can be repeated with similar success by substituting the generically or specifically 
described compositions, reactants and/or operating conditions of this invention for those used in the preceding exam- 
ples. 

[0112] From the foregoing description, one skilled in the art can easily ascertain the essential characteristics of this 
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invention and, without departing from the spirit and scope thereof, can make various changes and modifications of the 
invention to adapt it to various usages and conditions. 

[0113] While the invention has been described in connection with the preferred embodiment, it should be understood 
readily that the present invention is not limited to the disclosed embodiment Rather, the present invention is intended 
to cover various equivalent arrangements and is only limited by the claims which follow: 



Claims 

1 . A method of predicting a characteristic of a product to be prepared by a processing plant, comprising: 

(a) identifying an incoming material, said material having at least one measured property; 

(b) selecting at least one product to be prepared from said incoming material by said processing plant, said 
product having at least one desired characteristic; 

(c) selecting a process to be used by the processing plant to produce the product, said process having an 
effect upon the characteristic of the product; and 

(d) calculating a predicted value of the characteristic of the product using a predictive equation, 

wherein the predictive equation is generated by regressive analysis of an accumulation of data relating to 
(i) the or each measured property of the incoming material, (ii) the process and (iii) the characteristic of the product. 

2. A method according to claim 1 , wherein the incoming material is a petroleum feedstock, the product is diesel fuel, 
the process is admixing cetane improver, and the characteristic of the product is cetane number. 

3. A method according to claim 2, wherein the measured property is selected from hydrocarbon type, carbon content, 
hydrogen content, sulfur content, nitrogen content, API gravity, distillation range, cetane number, aniline point, 
heat content, kinematic viscosity, cloud point, pour point, and/or flash point. 

4. A method according to claim 2, wherein the measured property includes at least API gravity, T50 and at least two 
additional properties selected from T10, T90, aniline point, viscosity, cloud point, and aromatics content. 

5. A method according to claim 4, wherein either (a) the at least two additional properties are T90, aniline point, 
viscosity, and cloud point, (b) the at least two additional properties are T90, cloud point, and aromatics content, 
(c) the at least two additional properties are T10, T90, and cloud point, (d) the at least two additional properties 
are viscosity and aniline point or (e) the at least two additional properties are T1 0 and T90. 

6. A method according to claim 5, wherein the at least two additional properties are T90, aniline point, viscosity and 
cloud point and the predictive equation is a set of non-linear equations comprising: 

CNO = -0.07307 x T1 0 + 0.3080 x ANPT - 1 .1 52 x VISC40 + 0.4957 x D976 + 
0.1 836 x CLOUD +31 .57 



CN1000 = -0.09374 x T10 + 0.2738 x ANPT - 0.9512 x VISC40 + 0.5890 x D976 + 
0.2069 x CLOUD + 40.16 



CN2500 = -0.09681 x T10 + 0.2521 x ANPT - 0.9890 x VISC40 + 0.6700 x D976 + 
0.1943 x CLOUD + 42.38 



CN5000 = -0.09009 x T10 + 0.3084 x ANPT - 1.183 x VISC40 + 0.6461 x D976 + 
0.1695 X CLOUD + 41.30 
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CN7500 = -0.07602 x T10 + 0.3624 x ANPT - 1 .307 x VISC40 + 0.6495 x D976 + 
0.1371 X CLOUD + 35.77 

CN10000 = -0.06971 x T10 + 0.4190 x ANPT - 1.257 x VISC40 + 0.6033 x D976 + 
0.1086 x CLOUD + 33.75 

wherein 

CN is cetane number of the product at the noted concentration of cetane improver of from 0 to 10,000 ppmv, 
and intermediate values are interpolated between points, 
T10 is the temperature at which 10% of the product boils off, 
ANPT is the aniline point of the product, 

VISC40 is the viscosity of the product at 40°C according to ASTM D445, 
D976 is the cetane index of the product according to ASTM D976, and 
CLOUD is the cloud point of the product. 

A method according to claim 5, wherein the at least two additional properties are T90, cloud point and aromatics 
content and the predictive equation is a set of non-linear equations comprising: 

CN0 = 0.1 700 x T50 - 0.09022 x T90 - 238.4 x SPGR - 0.1 1 99 x FIAAROM + 
0.2187 x CLOUD + 237.8 

CN1000 = 0.1802 x T50 - 0.1062 x T90 - 256.1 x SPGR - 0.1 123 x FIAAROM + 
0.2419 x CLOUD 4- 259.8 

CN2500 = 0.1861 x T50 - 0.1073 x T90 - 264.7 x SPGR - 0.1213 x FIAAROM + 
0.2310 x CLOUD + 269.9 

CN5000 = 0.1928 x T50 - 0.1037 x T90 - 267.9 x SPGR - 0.1489 x FIAAROM + 
0.21 13 x CLOUD + 273.4 

CN7500 = 0.2079 x T50 - 0.09233 x T90 - 283.8 x SPGR - 0.1 700 x FIAAROM + 
0.1829 x CLOUD + 282.1 

CN10000 = 0.2163 x T50 - 0.08637 x T90 - 285.5 x SPGR - 0.1931 x FIAAROM + 
0.1574 x CLOUD + 281.7 

wherein 

CN is cetane number of the product at the noted concentration of cetane improver of from 0 to 1 0,000 ppmv, 
and intermediate values are interpolated between points, 
T50 is the temperature at which 50% of the product boils off, 
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T90 is the temperature at which 90% of the product boils off, 
SPGR is the speaUc (API) gravity of the product, 

FIAAROM is the a romatics content of the product according to ASTM D1319, and 
CLOUD is the cloud point of the product. 

5 

8. A method according to claim 5, wherein the at least two additional properties are T10 t T90 and cloud point and 
the predictive equation is a set of nonlinear equations comprising 

W CN0 = 0.05286 x T10 + 0.1329 x T50 - 0.07308 x T90 - 319.2 x SPGR + 0.1984 x 

CLOUD + 295.1 

, 5 CN1000 = 0.05529 x T10 + 0.1391 x T50 - 0.08732 x T90 - 332.5 X SPGR + 

0.2220 x CLOUD + 313.6 

20 CN2500 = 0.06523 x T1 0 4- 0.1357 x T50 - 0.08412 x T90 - 348.0 x SPGR + 

0.2085 x CLOUD + 328.0 

2 5 CN5000 = 0.091 78 x T1 0 + 0.1 1 82 x T50 - 0.06947 x T90 - 371 .4 x SPGR + 

0.1817 x CLOUD + 344.9 

30 CN7500 = 0.1203 x T10 + 0.1059 X T50 - 0.04557 x T90 - 403.8 X SPGR + 

0.1464 x CLOUD + 363.9 

35 CN 10000 = 0.1377 x T10 + 0.09925 x T50 - 0.03274 x T90 - 422.0 x SPGR + 

0.1158 x CLOUD + 374.5 



40 



wherein 



CN is cetane number of the product at the noted concentration of cetane improver of from 0 to 10,000 ppmv, 
and intermediate values are interpolated between points, 
T10 is the temperature at which 10% of the product boils off, 
T50 is the temperature at which 50% of the product boils off, 
4 5 T90 is the temperature at which 90% of the product boils off, 

SPGR is the specific gravity of the product, and 
CLOUD is the cloud point of the product. 

9. A method according to claim 5, wherein the at least two additional properties arc viscosity and aniline point and 
so the predictive equation is a set of non-linear equations comprising: 

CN0 = 0.281 1 x ANPT - 1 .030 x VISC40 + 0.651 9 x D976 

55 

CN1000 = 0.2403 x ANPT - 0.9091 x VISC40 + 0.7946 x D976 
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j 

CN250O = 0.2171 x ANPT - 1 .075 x VISC40 + 0.9147 x D976 



CN5000 = 0.2632 x ANPT - 1 .335 x VISC40 + 0.9365 x D976 



CN7500 = 0 3C48 x ANPT - 1 .322 x VISC40 + 0.9346 x D976 

10 

CN10000 = 0.3534 x ANPT - 1.259 x VISC40 + 0.9055 x D976 

wherein 

15 CN is cetane number of the product at the noted concentration of cetane improver of from 0 to 1 0 f 000 ppmv, 

and intermediate values are interpolated between points, 
ANPT is the aniline point of the product, 

VISC40 is the viscosity of the product at 40°C according to ASTM D445, and 
D976 is the cetane index of the product according to ASTM D976. 

20 

10. A method according to claim 5, wherein the at least two additional properties are T10 and T90 and the predictive 
equation is a set of non-linear equations comprising: 

25 CN0 = -0.5659 x T10 + 0.000458 x T50 x T90 - 0.000992 x (T90) 2 + 0.5261 x 

T10/SPGR + 0.4263 x T90/SPGR- 64.40 



30 CN1000 = -0.4433 x T10 + 0.000483 X T50 X T90 - 0.001215 x (T90) 2 + 0.4254 x 

T10/SPGR + 0.5364 x T90/SPGR - 81.53 



35 CN2500 = -0.4733 X T1 0 + 0.000468 x T50 x T90 - 0.001 226 x (T90) 2 + 0.4588 X 

T10/SPGR + 0.5420 x T90/SPGR- 79.54 



CN5000 = -0.5823 x T10 + 0.000399 x T50 x T90 - 0.001 130 x (T90) 2 + 0.5736 x 
T10/SPGR + 0.5086 x T90/SPGR- 73.52 



CN7500 = -0.6247 x T10 + 0.000359 x T50 x T90 - 0.001 125 x (T90) 2 + 0.6320 x 
T10/SPGR + 0.5245 x T90/SPGR- 79.91 



50 CN1 0000 = -0.7223 x T10 + 0.000346 x T50 x T90 - 0.001 046 x (T90) 2 +■ 0.7272 x 

T10/SPGR + 0.4876 x T90/SPGR- 74.34 

wherein: 

55 

CN is cetane number of the product at the noted concentration of cetane improver of from 0 to 10,000 ppmv, 
and intermediate values are interpolated between points, 
T10 is the temperature at which 10% of the product boils off, 
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T50 is Ihe temperature at which 50% of the product boils off, 
T90 is the temperature at which 90% of the product boils off, and 
SPGR is the specific gravity of the product. 

11. A method according to any one of the preceding claims, wherein the predictive equation is a set of non-linear 
equations derived by 

accumulating data relating the or oach said charactenstic of the product to (a) the or each measured property 
of the incoming material and (b) the effect of the process; and 

carrying out regressive analysis on the thus obtained data to correlate the said characteristic of the product 
to (a) and (b). 

12. A method as claimed in claim 1, wherein the product is diesel fuel, and at least one characteristic of the product 
is pour point. 

1 3. A method according to any one of the preceding claims, wherein step (d) is carried out on a computing device with 
appropriate software. 

14. A method according to any one of the preceding claims, wherein the cost of the incoming material and the market 
price of the product are known, and step (d) includes optimizing the profitability of the process. 

15. A method of manufacturing a product, comprising: 

(i) predicting a characteristic of the product by a method according to any one of the preceding claims; and 

(ii) manufacturing the product from the identified incoming material using the selected process. 

16. A method according to claim 15, comprising the further steps, repeated as necessary prior to step (ii), of 

changing the selected process and repeating step (i), to optimise the said characteristic. 

17. A computer program comprising computer program code means which, when executed on a computer, instruct 
the computer to carry out the steps of a method according to any one of the preceding claims. 

18. A computer readable medium having recorded thereon a computer program according to claim 1 7. 

19. A process for the value optimization of a plant which processes at least one incoming material and produces at 
least one product, comprising 

(a) assessing at least one property of the or each incoming material; 

(b) assessing the cost of the or each incoming material; 

(c) inputting said property and cost of the or each incoming material into a database; 

(d) determining at least one charactenstic of the or each product; 

(e) assessing the value of the or each product; 

(f) inputting said at least one characteristic and the value of the or each product into a database; and 

(g) calculating a value * optimized process for obtaining the or each product from the incoming material using 
a computing device accessing said database and calculating the cheapest way of producing the or each prod- 
uct, using a non-linear equation derived from regressive analysis of data from a history of properties of incoming 
materials, processes, and products. 

20. A process according to claim 19, wherein the plant comprises a petroleum products blending plant. 

21 . A process according to claim 1 9 or 20, wherein there are a plurality of properties of the or each incoming material 
which are assessed. 

22. A process according to claim 19, wherein there are a plurality of incoming materials and a plurality of products. 

23. A process according to any one of claims 1 9 to 22, wherein step (a) includes gathering information in at least one 
category selected from hydrocarbon type, carbon content, hydrogen content, sulfur content, nitrogen content, API 
gravity, distillation range, cetane number, aniline point, heat content, kinematic viscosity, cloud point, pour point, 
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and flash point. 

24. A process for the optimization of diesel fuel production, comprising 

(a) providing at least one database of diesel fuel stocks, additives and products, having a set of known prop- 
ertes; 

(b) providing a non-linear formula for the prediction of diesel fuel properties based upon a regressive analysis 
of the known properties collected from a series of samples; 

(c) providing a computing device connected to said at least one database; 

(d) providing computing instructions incorporating said formula for the prediction of diesel fuel properties; and 

(e) calculating the diesel fuel properties utilizing said computing device. 

25. A process according to claim 24, wherein the computing device is a computer and the instructions comprise com- 
puter software. 

26. A process according to claim 24 or 25, wherein the database is a spreadsheet of the set of known properties. 



27. A process according to any one of claims 24 to 26, wherein the database includes the price of the diesel fuel 
stocks, additives, and products and the cost of processing, and the software includes an optimizer, whereby the 

20 maximum profitability of the process may be calculated. 

28. A method of manufacturing a product, comprising: 

(a) optimizing the process parameters in a production plant by a process according to any one of claims 19 
25 to 27; and 

(b) manufacturing the product in the optimized production plant. 

29. A computer program comprising computer program code means which, when executed on a computer, instruct 
the computer to carry out the steps of a process according to any one of claims 19 to 27. 

30 

30. A computer readable medium having recorded thereon a computer program according to claim 29. 
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